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ABSTRACT Strong electron-electron interactions in graphene are expected to result in 
multiple-excitation generation by the absorption of a single photon. We show that the impact 
of carrier multiplication on photocurrent response is enhanced by very inefficient electron cooling, 
resulting in an abundance of hot carriers. The hot-carrier-mediated energy transport dominates the 
photoresponse and manifests itself in quantum efficiencies that can exceed unity, as well as in a char- 
acteristic dependence of the photocurrent on gate voltages. The pattern of multiple photocurrent 
sign changes as a function of gate voltage provides a fingerprint of hot-carrier-dominated transport 
and carrier multiplication. 



Graphene, a two-dimensional material with a gap- 
less electronic spectrum, possesses a unique combina- 
tion of optoelectronic characteristics [ijover a wide range 
of frequencies 0, Its high mobility enables high- 
speed photodetection^, [11 . Combined with optical trans- 
parency and gate-tunable carrier density (field effect), 
this makes graphene an attractive material for photonic 
and optoelectronic applications 0, One particularly 
appealing aspect of photoresponse in graphene, antic- 
ipated by theory but yet to be confirmed experimen- 
tally, is the possibility of multiple-excitation generation 
by a single photon|8[. This phenomenon, known as 
carrier multiplication (CM), is of fundamental impor- 
tance for optoelectronics, since many optoelectronic de- 
vices can achieve much higher efficiencies if operated in 
the CM regime. Although CM has been reported in 
nano-particles such as colloidal PbSe and CdSe quantum 
dots|9Mllj. at present the interpretation of these mea- 
surements remains controversial [l^ . 

Recent theory suggests that CM can be readily 
achieved in graphene [8| . Multiple-carrier production re- 
sults from impact ionization and Aug er- type processes 
induced by photoexcited carriers [13[ (see Figure [TJi). 
Fast carrier scattering, which dominates over electron- 



phonon scattering, is predicted by theorv [l4ll5l | and con- 



firmed by ultrafast dynamics observed in optical and ter- 
ahertz pump- probe studies [3, [13 ■ We distinguish the 
CM effect in graphene as 'intrinsic' from the effects in 
avalanche photodiodes, which operate at relatively high 
reverse bias, sometimes just below breakdown. 

Here we set out to understand the relation between CM 
and photocurrent response in graphene. We argue that 
the contribution of CM processes, which occur locally in 
the photoexcitation region, is enhanced by energy trans- 
port throughout the entire system area. In particular, 
because of slow electron-lattice relaxation which serves as 
a bottleneck process for electron cooling[l3|, photogener- 
ated carriers are thermally decoupled from the crystal 
lattice over length scales which, even at room tempera- 
ture, can be as large as ^ ~ 7//m. Thermoelectric cur- 
rents, arising in the presence of hot carriers, can lead to 
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FIG. 1: (a) Carrier multiplication in graphene via impact ion- 
ization of secondary electron-hole pairs, (b) Schematic of an 
optoelectronic device with two separately gated regions 1 and 
2, a laser excitation region positioned at the 1-2 interface, and 
a pair of contacts for collecting photocurrent. (c,d) Photore- 
sponse in homogeneous (c) and inhomogeneous (d) systems. 
Arrows indicate thermoelectric current due to photogenerated 
hot carriers. Maximal quantum efficiency is achieved for op- 
posite carrier polarities in regions 1 and 2, indicated by dif- 
ferent colors (d). 

a dramatic enhancement of photoresponse. 

As we shall see, these effects can have a direct impact 
on the quantum efficiency of photoresponse, namely the 
number of photogenerated carriers transmitted through 
the contacts per absorbed photon, q = A^ci/^ph- This 
quantity is a cumulative characteristic of the measured 
photoresponse, which depends on various effects occur- 
ring throughout the system, including the CM processes 
in the excitation region, as well as charge and energy 
transport from this region to the contacts. These pro- 
cesses are characterized by very different time scales: fast 
generation of carriers by the CM process is followed by a 
much longer "charge harvesting" stage dominated by the 
drift of carriers from the excitation region to the contacts. 

We emphasize that the hot-carrier regime arising under 
CM conditions is distinct from the effects of overall heat- 
ing. In the latter case, since the electron heat capacity 
is very small compared to the lattice heat capacity, only 
a small fraction of the absorbed photon energy, equipar- 
titioned between all degrees of freedom, would remain in 
the electron subsystem. As discussed in more detail be- 
low, this would result in a vanishingly small temperature 
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change, and suppression of the hot-carrier effects. 

In contrast, slow electron- lattice relaxation triggers 
thermal imbalance of the electron and lattice subsys- 
tems, amplifying the CM effects. The electron-lattice 
relaxation slows down for temperatures below the Debye 
temperature. Under these conditions, the extreme inef- 
ficiency of cooling mediated by acoustic phonons allows 
the carriers to remain hot during their entire lifetime be- 
fore reaching the contacts. Proliferation of hot carriers 
dramatically alters the nature of photoresponse. 

We exhibit the essential physics of photoresponse by 
considering a double-gated device comprising two regions 
with gate-tunable carrier densities (see Figure[T]D). As we 
will show, the quantum efficiency q for this system has a 
simple dependence on the local electrical conductivities 
cri.2 and chemical potentials ^1^2 in the two regions: 

where the factor a (0 < a < 1) describes the net frac- 
tion of the photon energy eo which is transferred to the 
electron system through photocxcitation and subsequent 
decay. Remarkably, under realistic conditions q does not 
depend on device dimensions and temperature, and can 
take values as high as g > 2 (sec detailed estimate below) . 

We characterize the electron system by an electron 
temperature, Td, which in general is different from the 
lattice temperature Tq. The electron energy distribu- 
tion is established via electron-electron scattering which 
occurs on a sub-picosecond time scale 1J| . Since these 
times are shorter than the electron-phonon timescales, 
the electron-lattice relaxation can be described by a two- 
temperature model. Crucially, the processes due to opti- 
cal phonons, which occur on relatively short times of sev- 
eral picoseconds 17[ , become quenched when the photo- 
generated carrier energies drop well below the Debye en- 
ergy, wd ~ 200 meV[19[. For the carriers with lower ener- 
gies, the dominant cooling process is mediated by acous- 
tic phonons, giving a slow cooling rate. The cooling rate 
mediated by the absorption and emission of an acoustic 
phonon was studied by Bistritzer and MacDonald (l8| . 
For cooling rate defined as dTc\/dt ~ — 7i(rci — To), 
Tel « To, they find 
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where = 20eV, p = 7.6 x lO^^g/cm^, and vp = 
lO^m/s are the deformation potential, mass density, 
and Fermi velocity for graphene. For realistic values 
/i = 100 meV and Tei ~ 300K, we arrive at timescales 
as long few nanoseconds. 

Slow cooling results in thermal decoupling of the elec- 
trons from the crystal lattice and energy transport me- 
diated by hot carriers. As we shall see, the effects due 



to energy transport dominate over the conventional pho- 
tovoltaic contribution to photoresponse. The hot-carrier 
mechanism and the photovoltaic mechanism of photore- 
sponse have very different experimental signatures. In 
particular, hot carriers manifest themselves in multiple 
sign reversals of photoresponse as carrier concentration 
is tuned by a monotonic sweep of gate voltage. Multi- 
ple sign changes do not occur in the absence of CM (see 
Fig 12]). Thus the pattern of photocurrent sign changes 
provides a fingerprint that can be used to experimentally 
identify the hot-carrier regime and the presence of strong 
CM. 

We describe the electric current in the hot-carrier 
regime through the local current density 



j = crE - er/n^(r)V[/g(r) + asVT^i 



(3) 



The first two terms describe the conventional photo- 
voltaic (PV) effect: primary photogenerated carriers are 
accelerated by the gate-induced electric field —VUg, and 
create a local photocurrent in the excitation region (here 
nx is the steady state density of photoexcited carriers and 
77 is the mobility at energy e ~ ■^hf). The redistribution 
of carriers associated with this local photocurrent sets 
up an electric field E = —\7{4> — M/e) that drives current 
outside the excitation region, reaching the contacts. The 
last term in Eq.© describes the contribution of energy 
transferred to electrons via multiple-carrier generation, 
which takes the form of a thermoelectric current driven 
by the electron temperature gradient. 

The quantities s and a in Eq.Q are the Seebeck coef- 
ficient and electrical conductivity, which depend on local 
carrier density and sign. The temperature profile can be 
found using the energy flux 



w= (0-^)j-nj-Kvrei, 
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where 11 = sT and k are the Peltier and the thermal con- 
ductivity coefficients. The values of ct, s and k depend 
on the microscopic scattering mechanisms. In the practi- 
cally interesting regime of disorder-dominated scattering, 
we have 
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where (...) denotes averaging over the energy distribution 
of carriers, and T(e) is the mean free scattering time for 
elastic collisions. 

We consider a simple model of photocurrent generation 
in graphene, based on Eqs.®,!]!]), which accounts for the 
multiscale character of photoresponse: fast carrier kinet- 
ics within a micron-size excitation region set up a pat- 
tern of local electron temperature and electric fields that 
drive current throughout the entire device. As shown in 
Figure [Ud, we consider a rectangular graphene sample of 
width W and length i, with a step in carrier density at 
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FIG. 2: Photocurrent map as a function of chemical poten- 
tials Hi and fj,2 for the device shown in Figure [T]3. Separately 
shown are the HC contribution (a) and the PV contribution 
(b), described by Eq.([7ll and Eq. (|lH) . respectively. Note mul- 
tiple polarity changes for the HC contribution as opposed to a 
single polarity change for the PV contribution. The scales of 
I(HC) and /(PV) have been calibrated to agree with the ratio 
calculated in Ea.(|17|l with L = 6/xm. Current slices along 
the dotted lines /ii = /x, shown for the HC contribution (c) 
and the PV contribution (d). The Seebeck coefficient s(^) is 
shown in (c). 

the interface between regions 1 and 2. Photocurrent is 
collected through two contacts placed at y = ±L/2. 

In this model, using the continuity relation V • j = 
and Eq.([3]), we express the photocurrent as 



/ / (.s(y)VT,i-a-i(y)e77n,V[/g)^ (6) 
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where R = ^ ^{y)dy is the total resistance, 

and the contacts are taken to be at equal potentials, 
I-L/2 ^y'^y ~ focus on the term s(j/)VTci in 

Eq.®, and for the time being ignore the other term. The 
latter contribution will be analyzed below and shown to 
be small. Approximating the dependence s{y) by a step 
that mimics the density profile, we express the hot-carrier 
(HC) contribution through the average increase of the 
electron temperature along the 1-2 interface 
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where si and S2 are the Seebeck coefficients in regions 1 
and 2, Tq is the lattice temperature, and the superscript 
'ave' stands for the value averaged over < x < W, 
y = 0. The spatial profile of Td must be determined 
from the heat transport equation 

V-W + 7Cei(Tei-To) =aeo7V, = hf, (8) 

where 7 is the electron-lattice cooling rate ([2]), Cei is the 
electron specific heat, and TV is the photon flux absorbed 



in the laser spot. Since typical spot sizes ^0.5 — 1 /im 
[2^ are smaller than other relevant length scales, such 
as the system size and cooling length (see below), the 
absorbed photon flux can be approximated as a delta- 
function source of hot electrons. 

Next, we analyze the HC photocurrent dependence on 
gate voltages of regions 1 and 2. Two simple cases are 
illustrated in Figures. [TJ;,d: the net HC current vanishes 
for a spatially uniform carrier density, and is maximized 
for a p-n interface. The full dependence on the chemical 
potentials //i,2, iflustrated in Figure [^t-: shows stronger 
photoresponse for fii and 112 of opposite sign. This is 
in agreement with the 'gate-activated photoresponse' ob- 
served in Ref. [20t in the presence of a gate-induced p-n 
junction. 

We model the dependence on the chemical poten- 
tials ^1^2, given by the factor ^^^^ in Eq.©, using 
R = ^^r^ and the Mott formulaic for the Seebeck 
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coefficient obtained from the non-interacting model ([S]), 
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where k^T max(/i, A). Here (j{fJ.) describes a linear 
dependence of conductivity on carrier concentration away 
from the Dirac point, with parameters Umin, the mini- 
mum conductivity, and A, the width of the neutrality re- 
gion. The dependence s{n), Eq.®, is in good agreement 
with the measurements of thermopower in graphene[22j. 

The dependence of photocurrent on /.(i and IJ.2 has a 
number of interesting features. Because of the depen- 
dence on s(/j,i) — 3(^2), the HC current ([7]) vanishes on 
the diagonal ^1 = ^2 and, in addition, on two hyperbo- 
lae /L<i/i2 = A^ that cross the diagonal. We shall refer 
to the latter as 'anomalous' polarity reversal. As illus- 
trated in Figure [5J;, this behavior can be traced to the 
non-monotonic character of the dependence s(^). In par- 
ticular, for any nonzero value oi ni excluding extrema of 
s(^), the dependence on /Z2 features two polarity changes. 
This is illustrated by the slice fii = /i* in Figure [2j:. At 
the nodes /i2 = m': the Seebeck coefficient satisfies 
s(/i*) = s(/i') = s(/i"), as shown by the horizontal dashed 
line in Figure [2j:. Hence the photocurrent has opposite 
polarity inside and outside the interval fi' < fi2 < fJ-" ■ 
Using a numerical evaluation of the integral in Eq.Q, 
we checked that multiple polarity reversal, 'proper' and 
'anomalous', as well as other qualitative features are in- 
sensitive to the photoexcitation spot size, surviving even 
for spatially uniform photoexcitation. 

A model similar to Eq.© was used in Ref.(23t to 
describe a laser-induced photo-thermoelectric effect ob- 
served in a heterogeneous system, a monolayer-bilayer 
interface. The measured photocurrent sign was consis- 
tent with Eq.([7]) but not with the PV effect (see below). 

The multiple polarity reversal in the dependence on fii 
and fi2 is unique for the HC mechanism. It is instructive 
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to make a comparison with the photocurrent response 
in the conventional photovoltaic (PV) regime where the 
primary photogenerated pair is the main contributor to 
photoresponse. This contribution is described by the sec- 
ond term of Eq.®, giving 
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a-^{r)er]n^{r)\7Ug{r)dxdy. (10) 



The integration simphfies when the size of the photocxci- 
tation spot is larger than the depletion length, Iq ^ Wd- 
Setting eUg{r) — /i(r), using the model dependence a{fi) 
from Eq.®, and replacing the integration over y by in- 
tegration over /z, gives 
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Result (|lip also shows photoresponse maximized in the 
presence of a p-n junction, i.e. for fii and fi2 of opposite 
sign. Thus the 'gate-activated photoresponse' of Ref.(20| 
cannot be used to distinguish the HC and PV contribu- 
tions to photocurrent. 

The difference between the two contributions to pho- 
tocurrent is most striking when /ii and fi2 have equal 
signs [see Figure [2^,b] . Since the polarity of the PV cur- 
rent is determined solely by the sign of field gradient 
VC/g, there is only one sign reversal occurring at fii = fi2- 
In contrast, the nonmonotonic character of s{fi) produces 
multiple polarity reversal of the HC contribution. Thus 
the polarity of photocurrent as a function of gate poten- 
tials offers a direct way to differentiate between the two 
mechanisms. 

To estimate the magnitude of the photocurrent, 
Eq. ([7]) , wc need to obtain the steady state profile of Tci 
from Eq.®. Since j has zero divergence, we can write 



V • W 
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(12) 



The first term, which is quadratic in j, can be ignored. 
The second term describes the Peltier coohng effect due 
to photocurrent passing through the 1-2 interface. Incor- 
porating it in Eq.® gives 

- V ■ (KVTei) + 7ai(Tei - To) = aeoN+3 ■ VH. (13) 

Since the spatial extent of the Peltier term j • VII is of 
order of the the p-n junction width, which is ^ 0.1 ^m in 
the state-of-the-art devices, it can be well approximated 
as a delta function source localized at the 1-2 interface. 

Eq. (jl3p can be conveniently analyzed using quantities 
averaged over the device width < a; < W, T'j5[™(j/) = 
^J^T,iix,y)dx, iV-™(y) =^ ^ Nix,y)dx, and 
transforming Eq. (|13p to a one-dimensional equation. For 
simplicity, we will consider the case when the laser spot 
is positioned on the 1-2 interface. Treating it as a delta 
function, we solve Eq. p^ by piecing together solutions of 



the homogeneous equation satisfying zero boundary con- 
dition at the contacts y = ±L/2, and performing match- 
ing of the boundary values at the 1-2 interface. 

Estimating the cooling length ^ = a/k/tC^ we find a 
large value ^ w 7 fim which exceeds L/2 for typical de- 
vice dimensions. Eq. p^ can thus be solved by approx- 
imating 7 « 0, yielding temperature profile ST^^°{y) = 
(1 — 2\y\/L)AT with the average temperature increase at 
the 1-2 interface 
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Since R cc L, AT grows linearly with system size, sat- 
urating when L/2 exceeds the cooling length ^. For 
not too high temperatures, fceT ^ max(//i^2, A), the 
second term in the denominator is smaller than the 
first term and can thus be ignored. Combining this re- 
sult with Eq.©, we evaluate the quantum efficiency as 
q = /(Hc)/(e/oW^-^yIo)- Using the Wiedemann-Franz 
(WF) relation e^K = ^k^Ta, and the Mott formula, 
Eq.®, we arrive at Eq.®. The result ® describes the 
realistic situation of large cooling length, £, ^ L. 

For a general system size L, the solution can be ob- 
tained as ST{y < 0) = Aisinh((y + \L)/£,i), ST{y > 
0) = A2 sinh((iL — y)/^2)- After matching the boundary 
values and derivatives at y = 0, wc obtain 
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This result agrees with Eq. p4)) for small system size 
L <C Ci,2- At large L, it describes saturation to the value 

AT ^ «eo?oiV-'o/('*i/ei + «2/6)- 

For an estimate of the numerical value of q, we use 
the factor a calculated in Ref.[^ which predicts aeo = 
M/cbTo with M = 4.3, Tq = 300 K, thus giving a « 0.07. 
Taking the chemical potentials at the minimum and max- 
imum of s(/i) described by the model ®, ^1^2 = ±A, we 
obtain q = MA:b'7o/2A. Taking a typical value for the 
neutrality region width for graphene on BN substrate, 
A < 300 K, we find q > 2. Thus, high quantum ef- 
ficiencies are feasible for realistic system sizes of up to 
5 — 10 /zm. 

Wc now proceed to estimate the relative strength of 
the HC and PV contributions to photoresponse. Using 
the WF relation and the Mott formula, Eq.®, we find 
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where we estimated the photoexcited carrier density as 
7ix = 2roA^ with tq the carrier lifetime. Near the Dirac 



point, estimating 
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for L ^ 5 — 10 /xm, where we used parameter val- 
ues: 7] — 10'^ cm^V~^s~^ , To = lpsl6|, the neutral- 
ity region width A « 100 meV estimated for graphene 
on SiO substrate 24 1 , a ~ 0.07 the photon energy 
eo ^ 1-5 eV, and the laser spot size Iq ~ I ^m psj . 

We therefore conclude that, due to very inefficient 
electron-lattice cooling (Eq.©) and efficient CM process, 
an abundance of hot carriers leads to a dominant HC con- 
tribution to the photoresponse. Furthermore, because 
the ratio /(hc) /^(pv) scales inversely with the square of 
the neutrality region width A^, we expect hot-carrier- 
related phenomena such as high quantum efficiency to 
become more pronounced for high mobility samples, e.g. 
graphene on a boron nitride substrate. 

To better understand the relation between the result 
PT)) and CM, it is instructive to consider the situation 
when energy is equally partitioned between all degrees of 
freedom, electrons and lattice (which would be the case 
for very fast electron-lattice relaxation). Crucially, the 
large difference between the electron and phonon specific 
heat values makes the lattice act as a nearly ideal heat 
sink. In this case, the value a « 0.07 used above, de- 
scribing the fraction of photon energy remaining in the 
electron subsystem, would be replaced by the heat ca- 
pacity ratio Cei/Cph, which is very small. In the ide- 
alized case of a sharp Dirac point and undopcd sys- 



tem, we have Cd/Cph 



IQ-"^. Although 



Cci is somewhat enhanced at finite doping, under re- 
alistic conditions it remains quite small. We estimate 
Cci/Cph « X il^/kBT), where keT < ii. Tak- 

ing ^ = 100 mcV, for temperatures lOK < T < 300K 
we obtain Cd/Cph ~ 10"^ — 10~^. This would reduce 
our estimate of the ratio /(hc)/-^(pv) by a large factor 
10 — lO'^, strongly suppressing ^(hc)- Hence, a dominant 
/(HC) signals that the electronic system and lattice arc 
out of equilibrium viz. the presence of CM. 

In summary, hot-carrier transport in the presence of 
photoexcitation leads to a novel type of photoresponse 
dominated by photo-thermoelectric effects. Multiple sign 
reversals of photocurrent over a monotonic sweep of gate 
voltage is a hallmark of this regime. The pattern of 
sign changes provides a clear fingerprint for experimen- 
tally identifying efficient carrier multiplication. Better 
understanding of nonequilibrium transport physics in 
graphene, and in particular of different pathways for cool- 
ing, is needed to develop a detailed picture of this new 
regime. This will open up broad vistas for both the ex- 
ploration of novel opto-electronic properties as well as 



the design of more efficient graphene photodetectors. 
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